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Abstract 

(N ■ 

^S) ' We examine the relative performance of algorithms for the calculation of curvature in spacetime. 

The classical coordinate component method is compared to two distinct versions of the Newman- Penrose 
7—i ■ tetrad approach for a variety of spacetimes, and distinct coordinates and tetrads for a given spacetime. 

^ ' Within the system GRTensorll, we find that there is no single preferred approach on the basis of speed. 

Rather, we find that the fastest algorithm is the one that minimizes the amount of time spent on 
' simplification. This means that arguments concerning the theoretical superiority of an algorithm need 

not translate into superior performance when applied to a specific spacetime calculation. In all cases it 
is the global simplification strategy which is of paramount importance. An appropriate simplification 
' strategy can change an untractable problem into one which can be solved essentially instantaneously. 

o> '■ 

1 Introduction 

I ' 

^ ' It is fair to say that tetrad calculations are generally considered superior to classical coordinate methods 
Zr', for the calculation of curvature in spacetime. Experiments by Campbell and Wainwright [Q now dating 
' back many years showed that tetrad methods are faster than coordinate methods by factors of 2 ^ 4. 
k>( I Even larger factors have been obtained by MacCallum |^ in a Euclidean context. Within the well known 
j_j ■ system SHEEP, for example, the advice to the beginner is to always use frame versions of the metric 
03 I (e.g. MacCallum and Skea |^ p. 23). On the commercial side, within the system MACSYMA2 Q the 
demonstration CTENS0R4 begins with an explanation that "frame fields" (orthonormal bases) allow the 
computations to run much more quickly. The demonstration calculates the bases components of the Ricci 
tensor for the Kerr-Newman spacetime in Boyer-Lindquist coordinates, and is a good place to begin our 
discussion. 



o 
o 



In Table 1 we have reproduced this demonstration within the system GRTensorll ||5|[|6| running under 
MapleV Release 3 Q, and have included the calculation of the Weyl tensor. The Table demonstrates 
some interesting properties. The theoretical advantage of the frame approach is clearly demonstrated in 
the Boyer-Lindquist coordinates (Column BKN). However, under the elementary coordinate transformation 
u = a cos 9 this advantage fails to deliver superior performance (Column BKNU). The importance of strategic 
application of simplification at intermediate steps is illustrated in Column BKNS. For this test simplification 
of components has been carried out only after the components of the Ricci and Weyl tensors are calculated. 
It is worth noting that without some optimization in the simplification strategy (e.g. post-calculation 
simplification only) this calculation cannot be executed in MapleV on a 32 bit machine. 

Clearly one could dismiss the findings in Table 1 if the implementation of the bases algorithms in GRTen- 
sorll were particularly incfhcicnt compared to coordinate methods. We do not believe this to be the case |8|. 

* e-mail: pollneyOastro .queensu. ca, musgraveSastro.queensu. ca, kevinOastro . queensu. ca, and lakeSastro . queensu. ca 



1 



Calculation 


BKN 


BKNU 


BKNS 


Riam 


6.7 


5.6 


30.0 


C(a)mc)(d) 


1.4 


0.8 


9.5 


Total 


8.1 


6.4 


39.5 


Rab 


8.4 


2.1 


18.9 


Cabcd 


21.1 


4.3 


67.5 


Total 


29.5 


6.4 


86.4 



Table 1: Average CPU time"'^ in seconds for the calculation and simplification of the bases com- 
ponents of the Ricci and Weyl tensors {R(a)(b) and C i^a){b){c)(d)) compared to the same for the 
coordinate components {Rab and Cabcd)- BKN refers to the Kerr-Newman spacetime in Boyer- 
Lindquist coordinates, and BKNU the same but with the transformation u = acos6. For both 
BKN and BKNU the simplification procedures have been structured for optimum performance^ '|^. 
For the Column BKNS the tetrad components of BKN are used, but simplification procedures are 
not applied to intermediate calculation steps, only to the final results. 



Rather, as we attempt to show in what follows, we believe that Table 1 reflects the fact that bases methods 
are not fundamentally superior to classical coordinate methods. We find that the most important criterion 
for speed is the minimization of the time spent on simplification. This underlines the importance of the 
choice of coordinates or tetrad in a computer algebra calculation and, more importantly, points out the fact 
that the user must be able to select the style of simplification which is most appropriate for the particular 
problem. 



2 Protocol for comparisons 
2.1 Choice of algorithms 

Within the framework of tetrad methods, the formalism of Newman and Penrose |^ has proven most useful 
for calculations in spacetime (see e.g. [0). Some of the earliest applications of computer algebra to relativity 
stressed the efficiency of this formalism (e.g. ||l|). McLenaghan ||l^ (see also Allen et al. |ll[) has emphasized 
two distinct approaches within this formalism. These are distinguished as the methods of Cartan and Debever 
in [l0|, and as Methods A and B in ||ll|, a notation which we adopt here. The methods are outlined in [l0[| 
and ]ll| with references and we do not repeat this material here|^. We simply note that Method A uses the 
definitions of Newman and Penrose explicitly, while Method B essentially uses definitions constructed so as 
to avoid inversion of coordinate indicies. 

In this paper we compare these two approaches to classical coordinate methods (suitably optimized). 



2.2 Basis for comparison 

The null tetrad formalism is sufficiently distinct from the classical coordinate approach that a basis for the 
comparison of the two methods is not clearly defined. In the classical approach the "curvature" of a spacetime 

^AU times_are in seconds as returned by the Maple V status function and are the average of four runs on a Sun Sparc 5 
(see Section 2.4). The maximum deviation from the average is less than 5% for times exceeding 2 seconds and about 10% for 
shorter times. 

■^We consider a worksheet (a sequence of calculation and simplification procedures) to be optimized when the execution time 
has reached a minimum. 

•^Due to their length, the complete text of worksheets used for these tests have not been included in this report (except for 
an example in Appendix ^) , however they have been made publicly available [ p^ . 

* The curvature component is consistently incorrect in pd] and [ |ll| . In particular, the coefficients of /xt and ua are -1 
and 1 respectively, not -2 and +2. This error is also present in the debever package in MapleV Releases 2 and 3. 
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is usually considered evaluated when the coordinate components of the Ricci and Weyl tensors have been 
evaluated. In the Newman-Penrose (NP) formalism it is the tetrad components of these tensors (the $s and 
4's) that display the "curvature" . The complication that arises in a comparison of such different methods is 
the fact that the natural output of each method is distinct. Now given the coordinate components, and the 
null tetrad, the tetrad components follow in the usual way One could then form a basis for comparison 
by defining the "curvature" as the tetrad components of the Ricci and Weyl tensors . This is the comparison 
used in Naturally, this puts the classical component method at a disadvantage since the extra sums 
involved are not a natural part of the method. In this paper we have tried to cover all possibilities by 
having both NP approaches output the tetrad components of the Ricci and Weyl tensors, and the coordinate 
approach output both the coordinate and tetrad components of these tensors. 



2.3 Choice of spacetimes 

To compare the performance of algorithms for spacetime calculations it is essential that a variety of spacetimes 
be considered, and that within a given spacetime different tetrads (coordinates) be examined. Campbell 
and Wainwright chose to examine the spacetimes of Griffiths Lcwis-Papapetrou [^, Bondi [ p^ , 
and Debever ||T^. More recently, the examination by Allen et al. (which is a comparison of the 

NP approaches) included these spacetimes (with a more general form of the Debever metric, the Debever- 
McLenaghan-Tariq metric p5[|) and also the plane wave, 2x2 decomposable, and static spherically symmetric 
spacetimes [|6| . We have found these last three spacetimes to be too simple since the associated calculation 
times are too short to form a reliable basis for comparison. They also include one form of the Kerr-Newman 
metric^ and a general tetrad (their Case 9). Here we examine the Kerr-Newman metric in a variety of forms. 
We do not include the tetrad 9 given in since it does not conform to the requirements of a null tetrad 
in the Newman-Penrose formalism iflql . 



2.4 Method of comparison 

The comparisons were made by way of 26 re-executable Maple worksheets, and the calculations were per- 
formed with GRTcnsorll j|] under MapleV Release 3 |0] (in the X- Windows interface) with patchlevel 3 on 
a Sun SPARC5 running SunOS 4.1.4 and equipped with a 75MHz CPU and 64Mb of RAM^ 

The Maple worksheets and associated input files used in these tests are summarized in Appendix The 
associated tetrads, and line elements are shown in Appendix The reproduction of tetrads is prone to 
errors, and references ||l|, |jl^ and ijlll all contain misprints in the tetrad components. Appendix |^ has 
been produced directly from the input files and so an error constitutes not a misprint, but an actual error 
in the input which would invalidate our conclusions in the case concerned. The worksheets are available for 
detailed examination and execution ||l7|| . 

The worksheets are constructed in the following way. Either the contravariant or covariant tetrad is 
loaded into GRTensorll from an input file. The choice determines NP algorithm to be used for the tetrad 
part of the calculation. The tetrad components of the Ricci and Weyl tensors are then evaluated, usually in 
simplified form^ (see Section The metric is then generated from the tetrad, and the covariant coordinate 
components of the Ricci and Weyl tensors are evaluated and simplified when necessary (again see Section ||). 
We have been careful to ensure that the metric, though generated from the tetrad, is presented in optimal 
form. Assuming that the tetrads are given, the worksheets then evaluate the tetrad components from the 
covariant coordinate components. These are then simplified to the exact form of the tetrad calculation. 

® Campbell and Wainwright |l| also examined the Kerr-Newman metric. However, their description involves intermediate 
hand simplification in this case and does not include comparison to the component method. In this Section we take the position 
that a fair comparison of computer algorithm efficiency can only be carried out if no hand simplification of tensor components 
is permitted. 

^This configuration was chosen for its reliability and reproducibility of CPU times, not its speed. By way of a comparison, 
a Pentium 133 is about 1.5 times faster, however the DOS/Windows implementation of MapleV Release 3 reports only integer 
CPU times. 

tensor component is considered to be fully simplified when its size, as measured in Maple 'words', reaches a minimum. 
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We believe that each worksheet has been fuUy optimized^. That is, the procedures (e.g. precalculation 
of the spin coefficients) and simplification procedures (type and order of simplification) at each step of the 
calculation have been constructed so as to present each approach in its best performance mode. It should 
be pointed out that partial optimization of a worksheet is straightforward (see Appendix ^), but the full 
optimization of a worksheet is a somewhat involved task. Full optimization is necessary if a real comparison 
of approaches is to be given. 

A sample worksheet is given in Appendix ^ along with the associated input file. In GRTensorll the input 
files contain no information beyond the components of the basis or metric. In particular, no simplification 
information is contained in the input file. The simplification strategy used is read from the worksheet. 



3 Comparisons 

Summarized in Table 2 are the results of our tests^ and the results available from the previous tests in 1^ 
and ||Tl[| . Before we discuss the comparisons it is appropriate to emphasize the importance of simplification 
procedures. The fastest calculation in the Kerr-Newman metric is the covariant NP tetrad calculation 
(Method B) in Row 9. If even only part of the simplification strategy is altered (e.g. the background 
simplification procedure used before the components are more fully simplified) the execution time can increase 
by a factor well over two orders of magnitude. The global simplification strategy is of paramount importance. 
A number of interesting points emerge from Table 2. 

i) It is appropriate to begin by comparing our results with previously published tests. Starting with the 
work of Campbell and Wainwright [Q, although the exact form of their output is unavailable, we observe 
a notable agreement for the ratio A/C and their results for the Lewis-Papapetrou and Bondi metrics as 
shown in Table 2. For the Griffiths metric our component times are somewhat fastei]^. Quite naturally, a 
concern at the time was the storage requirements for the calculations. They report storage requirements for 
the component method a factor 2^5 times that of the tetrad approach. Whereas storage is no longer the 
concern that it once was, we note that we have observed storage requirements for the component calculations 
only about 1.1 ~ 1.5 times that of the tetrad method. 

The paper by Allen et al. is concerned with the ratio A/D, that is, the relative performance of 
the two null tetrad methods. The exact form of their output is unavailable. In general we find that the 
performance of the Debever approach (Method B) is overestimated in ||ll|0. Although the central thesis of 
[pl| is the superiority of the Debever approach, this rests principally on their analysis of the Kerr-Newman 
metric. It is clear from Table 2 that this superior performance is tetrad (coordinate) dependent. 

ii) For metrics of a general type (Rows 3 through 8 in Table 2) we find that the standard NP approach is 
faster than the alternative proposed in and [pl]| . We find that this superior performance is not uniformly 
maintained if the general functions are replaced by specific ones. We have used the Kerr-Newman metric 
as an example, and as can be seen from Rows 9 through 14 of Table 2, the relative performance of the two 
tetrad approaches is highly dependent on the tetrads (coordinates). However, whereas the _De6e?;er approach 
can significantly outpace the NP approach (Rows 9 and 12), it is never far behind. 

iii) For metrics of a general type (Rows 3 through 8 in Table 2) both tetrad approaches are superior to a 
calculation of the tetrad components from the coordinate components. This is exactly as one would expect. 
The extra sums involved slow the coordinate approach down (compare Columns B and C as well as E and 
F ). Again, however, this superior performance is not uniformly maintained if the general functions are 
replaced by specific ones. Columns A/C and D/F for Rows 9 through 14 show that the classical approach 
can rival the tetrad approaches even for a calculation of the tetrad components. 

iv) It could be argued that the natural output of the coordinate approach, as regards the calculation of 
curvature, is simply the coordinate components of the Ricci and Weyl tensors. For metrics of a general type 
the standard NP approach retains its superiority over the coordinate calculation (Column A/B). This is not 

*The fourth component of for this metric given in jlj is wrong. We believe this to be a misprint which would not alter 
the time reported. 

^The errors in <I>i2 within the MapleV debever package affects only their result for the Debever-McLenaghan-Tariq metric 
(See footnote 4). 
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1. 


Coord's 


A 


B 


C 


D 


E 


F 


A/D 




A/B 


A/C 


|1 




D/E 


D/F 


2. 


Grif 


2.0 


1.9 


2.9 


1.6 


2.3 


3.5 


1.25 


1.65 


1.05 


0.69 


0.39 


0.70 


0.46 


3. 


L-P 


1.9 


2.8 


4.9 


2.6 


2.9 


6.1 


0.73 


1.44 


0.68 


0.39 


0.40 


0.90 


0.43 


4. 


Bondil 


3.5 


6.2 


10.5 


4.4 


6.2 


10.4 


0.80 


0.83 


0.56 


0.33 


0.27 


0.71 


0.42 


5. 


Bondi2 


1.5 


1.8 


4.7 


2.6 


1.8 


4.6 


0.58 


? 


0.83 


0.32 


? 


1.44 


0.57 


6. 


Deb 


11.1 


13.3 


39.8 


24.0 


13.6 


39.7 


0.46 




0.84 


0.28 






1.76 


0.60 


7. 


DMTl 


10.9 


13.5 


40.6 


22.3 


13.5 


40.7 


0.49 




0.81 


0.27 






1.65 


0.55 


8. 


DMT2 


35.3 


32.1 


109.0 


55.0 


32.6 


105.0 


0.64 


0.93 


1.10 


0.32 






1.69 


0.52 


9. 


KN-Eucl 


22.7 


13.1 


19.9 


6.3 


13.7 


20.0 


3.60 


5.80 


1.68 


1.14 






0.46 


0.32 


10. 


KN-Euc2 


54.4 


13.2 


17.2 


24.4 


13.1 


16.2 


2.23 




4.12 


3.16 






1.86 


1.51 


11. 


KN-BLl 


26.6 


27.5 


38.4 


30.0 


30.0 


34.6 


0.89 




0.97 


0.69 






1.00 


0.87 


12. 


KN-EFl 


32.5 


22.7 


30.2 


7.8 


22.7 


25.4 


4.17 




1.43 


1.08 






0.34 


0.31 


13. 


KN-BL2 


12.6 


6.7 


14.9 


18.2 


9.1 


15.8 


0.69 




1.88 


0.84 






2.00 


1.15 


14. 


KN-EF2 


41.0 


22.9 


38.9 


20.8 


22.8 


26.0 


1.97 




1.79 


1.05 






0.91 


0.80 



Table 2: CPU times and comparisons for optimized calculations. Columns A through F give CPU 
times in seconds^. 



Column 
A 

B 

C 

D 

E 

F 



total CPU time to generate and simplify the curvature components (<i>s and ^'s) 

from a contravariant tetrad using the standard NP approach ('Method A'), 
time for the simplified (covariant) coordinate components of the Ricci and Weyl 

tensors for the metric generated from the tetrad A. 
time (including B) to generate and simplify the tetrad curvature components from 

the coordinate components calculated in B, given the tetrad A. 
the time to generate and simplify the curvature components from a covariant 

tetrad using the modified NP approach of jl^ ('Method B'). 
time for the simplified (covariant) coordinate components of the Ricci and Weyl 

tensors for the metric generated from the covariant tetrad D. 
time (including E) to generate and simplify the tetrad curvature components from 

the coordinate components calculated in E given the tetrad D. 



The differences in Columns B and E are due to differences in the exact form of the metric 
generated from the tetrad. The spacetimes and form of the output are distinguished by the 
following abbreviations: 

Grif Griffiths metric [Q. Output in factored form. 

L-P Lewis-Papapetrou metric [|l^ . Output in factored form. 

Bondil Bondi metric Output in factored form. 
Bondi2 Bondi metric. Output in expanded form. 

Deb Debever metric |l^. Output in normal form (does not simplify further). 

DMTl Debever-McLenaghan-Tariq metric [|l^. Output in normal form (does not 
simplify further). 

DMT2 Debever-McLenaghan-Tariq metric in general form |15 . Output in normal 

form (does not simplify further). 
KN-Eucl Kerr-Newman metric |0. Output in factored form. 
KN-Euc2 Modified form of 9. Output in factored form. 

KN-BLl Kerr-Newman metric in Boyer-Lindquist coordinates |Q. Output in 
factored form. 

KN-EFl Kerr-Newman metric in advanced Eddington-Finkelstein coordinates [p^ . 

Output in factored form. 
KN-BL2 Kerr-Newman metric in modified Boyer-Lindquist coordinates using 

u — a cos 9. Output in factored form. 
KN-EF2 Modified form of 12. Output in factored form. 
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Calculation 


Mix 


Mixl 


Mix2 


Mix3 


R{a){b) 
C(a)(b)(c)(d) 


8.7 
8.5 


7.5 
1.1 


4.2 
7.4 


2.5 
1.1 


Total 


17.2 


8.6 


11.6 


3.6 


Rab 
Cabcd 


8.5 
52.5 


8.5 
13.7 


8.5 
52.5 


8.5 
13.7 


Total 


61.0 


22.2 


61.0 


22.2 



Table 3: Average CPU time in seconds for the calculation and simplification of the bases com- 
ponents of the Ricci and Weyl tensors {R(a)(b) and C(^a)(b){c)(d)) compared to the same for the 
coordinate components {Rab and Cabcd)- Mix refers to the standard 1-forms with trigonometric 
functions pO[ | and a time dependent basis inner product. The same inner product is used for 
Mixl but the trigonometric functions have been transformed away. For Mix2 a constant basis 
inner product has been used, and in Mix3 the trigonometric functions have been transformed 
away. In all cases the simplification procedures have been structured for optimum performance. 

true for the modified tetrad approach (Column D/E). Columns A/B and D/E indicate that at least in the 
Kerr-Newman metric the calculation of the coordinate components of the Ricci and Weyl tensors is usually 
faster than the calculation of the tetrad components. 

4 Discussion 

Our central conclusion is that the best algorithm, as regards speed, for the computer algebra calculation 
of curvature in spacetime is the one that minimizes the amount of time spent on simplification. This 
underlines the importance of the careful choice of coordinates or tetrad in a computer calculation and, more 
importantly, demonstrates that the user must be able to style the global simplification strategy in a manner 
most appropriate for the particular problem being studied. An appropriate simplification strategy can change 
an untractable problem into one which can be solved essentially instantaneously]^. Our comparisons (Table 
2) indicate that there is no uniformly superior algorithm. In the development of these comparisons we 
observed that the differences between procedures optimized with respect to the global simplification strategy 
for the procedure were less than the variations within a given procedure for different simplification strategies. 

Although we have not found any algorithm to be uniformly superior, there certainly are cases where the 
appropriate algorithm stands out. A case in point is the mixmaster spacetime |l9||^ Here an appropriate 
choice of basis removes any angular dependence in the bases components of both the Ricci and Weyl tensors. 
The tetrad approach (not a null tetrad in this case) would certainly be expected to outperform a coordinate 
calculation in this case. This is confirmed in Table 3. We have considered both a constant basis inner product 
(for Mixl and Mix3) and a time dependent one (for Mix and Mix2). The coordinate transformations are 
simply Q — cos 9 and 4" = sin-0. It is clear that the classical coordinate calculation is no match for the basis 
approach exactly as one would guess. Interestingly, it is the Weyl tensor calculation that improves under 
coordinate transformation, and the Ricci tensor under change in the basis inner product. 
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A Simplification strategies for GRTensorll in MapleV 



We have found that the global simplification strategy is of central importance for the computer algebra 
calculation of curvature in spacetime. Summarized here are some general rules of simplification appropri- 
ate to GRTensorll running under MapleV. Although the general philosophy described here is more widely 
applicable, our comments are specific to the system used and so we have relegated this information to an 
Appendix. 

(i) There should be some default simplification procedure applied to every step of a calculation. Failure to 
do this can make a simple calculation intractable. The Maple routine normal is a good starting point 
for the default. If the calculation involves exponentials (e.g. the Bondi metric) the routine expand 
may be more appropriate. Only if very general functions are involved is it appropriate to consider no 
default simplification. In general cases this may be the optimal choice. 

(ii) For both coordinate and tetrad calculations the removal of trigonometric and like functions via ele- 
mentary coordinate transformations will often improve performance. 

(iii) Simplification of the metric tensor or tetrad components before further calculation will improve per- 
formance. 

(iv) Precalculation and further simplification of the spin coefficients (and their complex conjugates) will 
improve performance only in more complicated cases. The same holds for the Christoffel symbols in 
the coordinate approach. 

(v) For further simplification after an object has been calculated, the Maple routine simplify is seldom a 
good first choice. The routine expand followed by factor is often more appropriate. If the situation is 
suHiciently general (e.g. the Debever-McLenaghan-Tariq metric) there will be no further simplification 
if normal has been used as default. 

(vi) If complicated functions are involved, it can be advantageous to substitute the explicit forms of the 
functions after a more general calculation is completed. 

(vii) When a calculation is proceeding slowly, it should be halted, the simplification strategy altered, and 
the worksheet reexecuted. 

For most situations these general rules will give adequate performance, and reduce the calculation of 
curvature for even complex spacetimes to an essentially trivial exercise. Usually, it is the answer that is of 
interest and not the fact that the simplification strategy is optimal. When optimal simplification strategies 
are the prime concern ( as in this paper ) the problem is more involved because of the large number of 
simplification procedures available^ and the size of the resultant parameter space to be explored. 

^■^ GRTensorll provides a menu of 12 distinct commonly used predefined simpUfication routines with the ability to introduce 
customized constraints and simplification routines. Any single parameter routine can be applied with gralter, and any multiple 
parameter routine can be applied with the command grmap. 
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B Worksheets and input files 



This Appendix is simply a bookkeeping operation which relates the input files necessary to run the 
worksheets. The first three worksheets are associated with Table 1, the last four with Table 3, and all 
the rest with Table 2. The worksheet appc.ms generates Appendix ^ from the input files. Abbreviations 
which appear twice in the last column of this Table refer to contravariant and covariant tetrad components, 
respectively. 

Worksheets used to generate Table 1: 



Test 


Worksheet 


Input file 


BKN 

BKNU 

BKNS 


bknl .ms 
bkn2 . ms 
bknS . ms 


bkn . mpl 
bknu . mpl 
bkn . mpl 



Worksheets used to generate Table 2: 



Test 


Columns A-C 


Columns D-F 


Worksheet 


Input file 


Worksheet 


Input file 


Grif 


grif f 1 .ms 


npupgrif .mpl 


grif f 2 .ms 


npdngrif .mpl 


L-P 


lewisl .ms 


npuplew.mpl 


lewis2 .ms 


npdnlew.mpl 


Bondil 


bondil .ms 


npupbon . mpl 


bondi2 .ms 


npdnbon . mpl 


Bondi2 


bondil .ms 


npupbon . mpl 


bondi2 .ms 


npdnbon . mpl 


Deb 


debl .ms 


npupdeb . mpl 


deb2 .ms 


npdndeb . mpl 


DMTl 


dmtsl .ms 


npupsdmt .mpl 


dmts2 .ms 


npdnsdmt . mpl 


DMT2 


dmtl .ms 


npupdmt . mpl 


dmt2 .ms 


npdndmt . mpl 


KN-Eucl 


knal .ms 


npupknl .mpl 


kna2 . ms 


npdnknl .mpl 


KN-Euc2 


knbl .ms 


npupkn2 . mpl 


knb2 . ms 


npdnkn2 . mpl 


KN-BLl 


kncl .ms 


npupknS . mpl 


knc2 .ms 


npdnknS . mpl 


KN-EFl 


kndl .ms 


npupkn4 . mpl 


knd2 . ms 


npdnkn4 . mpl 


KN-BL2 


knel .ms 


npupknS . mpl 


kne2 .ms 


npdnknS . mpl 


KN-EF2 


knf 1 . ms 


npupkn6 . mpl 


knf 2 .ms 


npdnkn6 . mpl 



Worksheets used to generate Table 3: 



Test 


Worksheet 


Input file 


Mix 
Mixl 
Mix2 
Mix3 


mixmr .ms 
mixmrl .ms 
mixmr2 .ms 
mixmrS.ms 


mix. mpl 
mixl .mpl 
mix2 .mpl 
mixS.mpl 



9 



C Tetrads and metrics produced from the input files 



The following sections list tetrads and metrics used as inputs for the tests listed in Tables 1-3. This list has 
been produced directly from the input files which were used in the tests (listed in Appendix ^) and converted 
to using MapleV's I^TJ^X output facility with only minor modifications to improve readability. 

C.l Kerr-Newman "frame field" (Table 1) 

For this set of tests (whose output is given in Table 1), the Kerr-Newman spacetime is described by a 
frame consisting of four independent covariant vector fields whose inner product is the constant matrix 
'7(a)(b) — diag(— 1, —1, —1, 1). The basis vectors and corresponding line element for each case are given below. 



Spacetime: Kerr-Newman (Boyer-Lindquist coordinates) (BKN and BKNS) ||T8| 



Input file 



Frame 



bkn .mpl, 



W4a 



+ a2 cos(6')2 

0, +a2 cos( 61)2, 0,0 
( r2 + a2 ) sin( e ) 



,0,0,0 



0,0, 



0, 0, 



a sin( 9 ) 



^/r^'^~aFcx)s{W)^ ' y^r2~t~a2'cos{^72' 
a ■\/r2 — 2 m r + a2 + (52 gj]^(| q ^2 y'r^ ~ 2mr + a'^ + 



y/r"^ + a'^ cos( 9 )2 
Corresponding line element: 
2 (r2 -|-a2cos(6')2) dr^ 



ds 



y/r^ + a2 cos( 9 )2 

, , , ^, + (-r^ - cos{9f) d9^ +i-l + cos{9)) 

(cos(6l) + 1) (r2 a2cos(6i)2 -2a2TOrcos(6')2 + a4cos(6i)2 

+ a^Q'^C0B{9f + 2a^mr + r^-a^Q'^ + r^a^^ d 4>^ J {r^ + cos{9f) 

2a(-l + cos(6'))(cos(0) + l)(-2mr + g2) d(j) dt 

r2 + a2 cos( 9 Y 
(a2cos(0)2 + r2 -2mr + Q2) dt^ 

r2 +a2cos(6')2 

Spacetime: Kerr-Newman (Boyer-Lindquist coordinates, u = a cos 6) (BKNU) 



Input file 



Frame 



bknu . mpl 



W2a 



W3a 



W4a 



Vr2+t? 



•\/r2 — 2mr + a2 + i 



,0,0,0 



0, - 
0, 0, 
0,0,^ 



vr2+i? 
y/ — U' 



,0,0 



( r2 + ) yj cP- — V? \J d?- — v? 



\/r2 + - 



+ m2 a 

^r2 — 2 m r + a2 + (32 Q,2 _ y2 ^ ^j.2 _ 2 m r + a2 + Q2 



\/r2 + ' 



^7*2+' 
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Corresponding line element: 

2 ( + ) d (r^ + u'^) du^ 
as — ^ ^ — r — 7 + (a — u) (a + u) 

( —r'* ~ (P' ~ V? — 2 c? mr ^' 2mr V? — c? iP' ^ c? Q"^ — Q"^ V? ) d (j)^ /( 

2 2N, 2^ (a-u)(a + u)(-2mr + (32)c?0c?i ^ (u2 + r2-2mr + (5^)dt^ 



( + ) a 



C.2 Null tetrads (Table 2) 

This section lists the set of null tetrads for the test cases used to generate Table 2. For each spacetime, 
four forms of input were used. The times listed in Column A of Table 2 are obtained using a contravariant 
null tetrad as input. In Columns B and C the metric is calculated from this tetrad and used for subse- 
quent calculation. In Column D a covariant tetrad is loaded, and in E and F, its corresponding metric is 
used. Though the metric of column E is, of course, equivalent to that of _B, there are often differences in 
representation which can in principle alter calculation times (though in practice we have found this effect to 
be minimal). Thus the line elements used as inputs for each tests are listed along with the tetrad used to 
calculate them. 



Spacetime: Griffiths (Grif) P 



Input file 



Contravariant tetrad 



npupgrif .mpl ^ = [0,1, 0,0] 
1, 0, 0, 0] 



n 

m° 



1 (-27a 2a t; + ay), ^ e^ "^^"^ ( ) ) ^3 ( -/ + 2 a u + ay ), 

o o 

lg(-2/a(«+.))^^ lg(-2/a(n+.))^(j_2a^,„2at.) 

Ig(2/a(«+t,))y3(j_^2aw + ay), -e(2^''("+"))V3(/ + 2au + ay), 
3 3 

1 g(2/a(«+^,))y3^ie(2-f«(«+^'))73(„/„2au-2at;) 

Corresponding line element: 

ds =2 du d V ^ 2 a {y — 2v) du d x +2 du d y + 2a{2u — y) dv d x 
+ 2 dv d y + ^—4 a'^yu + 2a^ y^ — Aa^vu — Aa?vy — 6a^ — ^ — 6 dx^ 

— 2a{u + 2y + v) d x d y + — dy^ 



Input file 



Covariant tetrad 



npdngrif .mpl la — 

na = 



1, 0, a{2u-y), 1] 
0,l,a{2v-y),l] 

0, 0, V3e(-2^"("+'')) (l + 2/a(ii + w)),-i/\/3e(-2^"("+")) 

0, 0, V3e(2^'^(«+'')) (l-2/a(ii + t;)),i/%/3e(2J'^(«+'')) 
Corresponding line element: 

rfs^ — 2 du d V — 2 a {y — 2v) du d x +2 du d y + 2a{2u — y) dv d x 
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+ 2 dv d y + y—A a? y u + 2 y"^ — 4:0? v u — 4:0? v y — Q a? v? — - — Q a? ) dx 
— 2a{u + 2y + v) dx d y + - dy^ 



Spacetime: Lewis-Papapetrou (L-P) 



Input file 



Contravariant tetrad 



npuplew.mpl 1° 



,0,0,1 



s(x,y)^ 



V2 



V2 



2 r(a;,y)e(-*'(^'2')) 
2 r(a;,j/)e(-''(^'2')) 



e(k(3:,j/)-s(a;,a)) 

%/2 



0, 



Corresponding line element: 



,0,0 
,0 



ds^ — 



dt^ 



^^{x,y)e<^^'y) dt d z _(^(,(.,,)_3(.,,)))2 ^^^y 



(e(-s(^.a)) )2 e^-'^t^^'S^)) 
(g(-s(:E,a)) _^(^^y)gs(:r,y) ) ( e( ) ) !■( x, ?/ ) + w( a;, y ) e"( ) ) rfz^ 



Input file Covariant tetrad 



npdnlew .mpl la 



ie^(^'!') ^/2,0,0,- ir(a;,y)e(-"(^'2')) \/2 + i w(a;,?;)e"(^'2') ^2 

i e''( ) ^ ^ , , i r( y ) e( -'^^ ^'2' ) ) V2 + i w( a;, y ) e''^ ^'2' ) V2 

0, 0, - ie('^(=^'2^)-"(^'S')) ^/2,0 

0, - ie(k(=^.a)-s(a;,a)) ^^0,0 
Corresponding line element: 

ds^ ^ ei^<--v)) dt^ +2e(2-(^'^»w(x,y) dt d z „ e(2k(-,y)-2s(.,,)) ^ ^ 
_ ( e( ~< '^■y)^r{x,y)-w{x,y) e< ^ ) ( e( ^^''■y ) > r( x, y ) + w( x, y ) e"( ) ) d z^ 



Spacetime: Bondi (Bondil, Bondi2) [|T|] 



Input file 



Contravariant tetrad 



npupbon.mpl l" = [e( "Qt'^"'^) ' , 0, 0, 0] 

1 e(-Q('''"'®)) V(r,u,6l) 



^ g( -Q( r,«,fl ) ) ^ g( -Q( r,u,9 ) ) U( r, U, ) , 



1 \/2e(-'^(''^"''')) 1 /e^('''"'^) V2 
" 2 r ' " 2 rsin(^) 



0, 0, 



1 V2e(-'^('''"'^)) 1 /eT('''"'^) 



2 r ' 2 r sin( 6* ) 

Corresponding line element: 
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Input file 



Covariant tetrad 



npdnbon.mpl = 



0, eQ( 0, 0] 

' r ' ' 

0, - irU(r,it,6l)eT('''"^'') ^/2, ^ r e^('^'"'^) V2, ^ /r sin( 6* ) e^ "T^'''"^^ ^ ^ \/2 

0, - irU(r,u,6l)eT('''"'^) \/2, ^ re'^^''-"-^ ) V2, - ^ / rsin( 6* ) \/2 



Corresponding line element: 

ds^^2i eQ( ) )^ . . . . + f 2 (^''^'-''rnr,u,e) ^ ^, 



2 / 7(r,«,e) n2- 



vir,u,ey 

2 a\2 J J,2 



Spacetime: Debever (Deb) |[T5| 



Input file 



Contravariant tetrad 



npupdeb .mpl 1°- = 



1 



2 -L(x,y)P(a:,y)+M(x,y)N(x,2/)' 



2 -L{x,y)Fix,y) + Mix,y)Nix,y)' 2 S(x,2;) 

1 P{x,y)V2 

2 -~L{x,y)F{x,y) + M{x,y)N{x,y)' 

1 N(a;,;/)72 1 Y{y)V2 

2 -L(x,y)P(x,y) + M(x,y)N(x,y) '2 S(a:,y) ' 

1 M{x,y)V2 

2 -L{x,y)P(x,y) + M{x,y)N{x,y)' 
1 L(a;,y)V2 

' 2 -L{x,y)P{x,y) + M{x,y)N{x,y) 

1 M(x,y)72 

2 -L(a;,y)P(a;,y) + M(rE,y)N(a;,y)' 
1 L(a;,y)y2 



1 IX{x)V2' 



2 Six,y) 



lIX{x)V2' 



2 -Ux,y)P{x,y) + M{x,y)N{x,y)' '2 S(x,y) 
Corresponding line element: 

ds^ = {L{x,y)-N{x,y)){L{x,y) + N{x,y)) dt' 
+ 2{-P{x,y)N{x,y) + M{x,y)L{x,y)) dt d z 



iP{x,y)-Mix,y))iP{x,y)+Mix,y)) dz 



2 S{x,yf dy^ 



S{x,y)^ dx^ 



Input file 



Covariant tetrad 



npdndeb . mpl la — 



iL(x,y)y2,lM(x,y)V2,,i^%^,0 
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'-lN,..„v^,-lp,.,,,v^,0,ii5ii^ 

4N(..„V2,-lp,.„A0,-ii^^ 
Corresponding line element: 

ds^ = mx,y)-N{x,y))mx,y) + Nix,y)) dt^ 
+ 2i-Pix,y)N{x,y) + Mix,y)L{x,y)) dt d z 

+ {Mix,y) - P{x,y)) iP{x,y) + M{x,y)) d - 



S{x,yf dy^ S{x,y)^ dx'^ 



X(x)2 



Spacetime: Debever-McLenaghan-Tariq (DMTl) 
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Input file Contravariant tetrad 



npupsdmt .mpl 1°" = 



V{w,x)V2 



N(w,a;) \/2 



2 \.{w,x)V{w,x) -M{w,x)n(w,x)' 2 l.{w,x)V[w,x) - M{w,x)'H{w,x)' 
1 V2 



0, -- 



2 S{w^x] 



V{w,x)V^ 



N(w,a;) ^/2 







2 L(w,a;)P(w,a;) - M(w;,a;)N(u;,a;)' 2 L( tu, a; ) P( tu, a; ) - M( a; ) N( w, a; ) ' 
1 \/2 



2 8(^,3;' 
1 



M(w,a;) V2 



L(w,a:) %/2 



2 L(-u;,a;)P(w,a;) - M(w,a;)N(w,x) ' 2 L( w, a; ) P( x ) - M( u;, a; ) N( x ) ' 
1 /\/2 



2 R(u;,a;) 
1 



M(w,a;) V2 



1 



L(w,a;) V2 



,0 



2 L(u;,a;)P(u;,a;) - M(w,a;)N(w,x) ' 2 L( w, x ) P( w, a; ) - M( u;, a; ) N( a: ) ' 

1 

2 R(w,a;) 
Corresponding line element: 

rfs^ = ( L( w, a; ) - N( a; ) ) ( L( w, a: ) + N( w, a; ) ) d 
+ 2 (-P(u;,a;)N(w,a:) + M( a; ) L( w, a: ) ) du d v 

+ (M( w,a;) - P(w,a:) ) (M(u;,a;) + P(w,a:) ) dv^ -'R{w,xY dw^ -S{w,xy d x^ 



Input file 



Covariant tetrad 



npdnsdmt .mpl la 



i V2L(it;,x), i%/2M(u;,a;), 0,is(w,a;) V2 
i V2 L( a: ) , i x/2 M( a; ) , 0, - i S( a: ) 72 

- i \/2 N( u;, a; ) , - i \/2 P( u;, a; ) , i / R( w, a; ) V2 , 

- i \/2 N ( u- , a; ) , - i V2 P ( , a; ) , - i / R ( w , a; ) V2 , 
Corresponding line element: 

ds = ( L( w, a; ) — N( w, a; ) ) ( L( ui, a: ) + N( w, a; ) ) du"^ 



14 



+ 2(-P(w,x)N{w,x) + M(w,x)L{w,x)) du d v 

+ {M{w,x) - P{w,x)){M{w,x) +P{w,x)) dv^ -Yi{w.xY dw^ -S(w,a;)2 dx^ 



Spacetime: Debever-McLenaghan-Tariq (modified) (DMT2) ||T5| 



Input file 



Contravariant tetrad 



npupdmt . mpl 1°" — 



P( u, v,w,x) \pl 



m 



,0, 



1 V2 



1 

2 1j{u,v,w,x)P{u,v,w,x) — M{u,v,w, x)N{u,v,w, x] 
1 N{u,v,w,x)V2 

2L{u,v,w,x)P{u,v,w,x)—M{u,v,w,x)'N{u,v,w,x)' ' 2 S{w,x) 

1 P{u,v,w,x) ^/2 

2 L( u, v,w,x) P( u,v,w,x) — M( u, v,w,x) N( u,v,w,x)' 
1 N(u,w,w,a::)\/2 



,0,- 



V2 



2 1j{u,v,w,x)P{u,v,w,x) —M{u,v,w,x)'N{u,v,w,x) ' ' 2 S(w,a;) 

1 M(u,i;, \/2 

2 L( u, w, a; ) P( M, i), w, a; ) — M( u, w, w, a; ) N( u, w, a; ) ' 

1 L(M,u,w,a;) V2 1 I V2 ^ 

2 1j{u,v,w,x)P{u,v,w,x) — M( u,v,w,x)N{u,v,w,x) 2 R( w, a: ) ' 
1 M{u,v,w,x)\/2 



1 



2 L( u, V, w, a; ) P( M, u, w, a; ) — M( u, w, w, a: ) N( u, w, w, a; ) ' 

L(M,t;,w,a;) \/2 1 /\/2 



2 L( It, u, w, a; ) P( u,v,w,x) — M( u, w, a; ) N( u, w, w, a: ) '2 R( w, a; ) 
Corresponding line element: 

ds^ = — (N(u,w,w,a;) — L(u, w,i(;,a;))(N(M, t;,w,a;)+L(u,w,w,x)) du^ 
+ 2 ( — P( u, w, w, a; ) N( u, w, w, a: ) + M( u,v,w,x)'L{u,v,w,x)) du d v 
— {P{u,v,w,x) — M( u,v,w,x)) {P{u,v,w,x) +M{u,v,w,x)) dv^ 
— R{w,x)'^ dw^ ~S{w,x)^ dx^ 



,0 



Input file 



Covariant tetrad 



npdndmt . mpl la 



^ V2L( u,v,w,x) , - %/2M( u,v,w,x) ,0, -S{w,x) \/2 
^V2L{u,v,'w,x), i \/2M(u,w, w,a;) ,0, -is(w,a;)y2 

- i V2 N( u, w, a: ) , - i x/2 P( u, w, a; ) , i / R( a; ) V2 , 

- i \/2 N( u, w, w, a: ) , - i V2 P( u, w, w, a; ) , - i / R( w, a: ) \/2 , 
Corresponding line element: 

ds = — ( N( u, w, w, a; ) — L( u, w, w, a; ) ) ( N( M, u, a: ) + L( u, w, w, a; ) ) d'l 
+ 2 ( — P( u,v,w,x)N{u,v,w,x) + M( w, t;,w,a;)L(w, u,w,a;)) rfu d v 
+ {M{u,v,w, x) — P{u,v,w,x)) {P{u,v,w, x) + M{u,v,w,x) ) dv^ 
-R(w,a;)^ dw^ -S(w,a:)^ da:^ 



Spacetime: Kerr-Newman {x,y,z,t) (KN-Eucl) 



Input file 



Contravariant tetrad 
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npupknl.mpl P 



V2a 



2 V.t2 -2Mx + q2 + Q2 ^ j.2^a2y2 '2 + a? y"^ sj x"^ - 2 M x + a? + Q"^' 



2 „2 



,0 



(a;2+a2)v/2 



\/2a 



.2 v^a;2 _ 2 M X + a2 + Q2 ^ y"^ ' 2 ^a;^ + a2 y V^^^ - 2 M a; + + g2 ' 
1 \/2v'a;^-2Mx + a2 + Q2 



m" = 



■\/a;2 + t/^ 



^/2 



1 7\/2yr^ 



2 a/x^ + y2 
^ 1 /^/2v/^^ 
2 -J x^ ■\- o? 



2 Y^a;2 + a2 ' 2 ^ -\- y^ y/ \ — y^ 

1 V2aVl - 1 

2 + a2 2/2 ' 2 y'a;2 + a2 ?/2 _ y2 

Corresponding line element: 

_ (a^;/^ + a;2-2Ma; + Q2) ^ ^ g ( _2Ma; + Q2 ) ( -1 + y ) + 1 ) ^ 



a;^ + o"' 1/ 



2 1,2 



a;2 + o2 «2 

+( a;^ + x2 a2 + 2 a2 M a; - o2 Q2 + ^2 ^2 ^2 _ 2 ^2 ^ ^ ^4 ^2 ^ y2 q2 ^ 
[ i+y)[y+L)az /[x+ay) ^2 _ 2 M x + + ^ {-1 + y) {y + 1) 



Input file 



Covariant tetrad 



npdnknl . mpl la 



1 72 v'a;^ - 2Ma; + a2 +Q2 _ l y2a v'a;2 - 2Ma; + a2 + Q2 ( 1 _ ^2 ^ 

'~ 2 



"a = 



2 y/x"^ +0? 2/2 

1 \/2 \/x2 + a2 y2 

2 ^x2 -2Mx + a2 + g2 
1 V^A/a;2-2Ma; + a2 + Q2 



-\/a;2 + a2 ?/2 



,0 



\/2^ 



-v/a;2 +a2 ^2 



1 V2 g v^a:2 _ 2 M x + a2 + g2 ( 1 - y2 ) 
2 



2 ,,2 



a;^ + 2/ 



,0 



mo = 



-^a;2 + 0? 2/2 

1 /v^V^2-:r^2^ 



2 v/a;2 - 2Ai".T + a2 + Q2 
'1 V2a^^g 1 V2VT^(a;2+a2) 
2 •\/!r2"+"a2y2 2 ^Jx'^ + a? y^ 

lV2a^/T^f_ 1 \/2VT^(a;2 + a2) ^ 1 iV2V^fTaV 

2 v'a;^ + a2 2/2 ' 2 ^a;2 + 2/2 ' '2 



'0'-2 



Corresponding line element: 

, 2 (a22/2 +x2 -2Ma; + Q2) dt^ 
ds = '■ — 



a;2 + a2 y2 



,(-2Mx + Q2)(_i + y)(y + i) dt d z 
x^ + 2/2 



+ {x'^ + x^ + 2a^ M X ~ + y^ x^ - 2a^ y^ M x + 2/ + a y Q ) 



(-l + 2/)(2/+l) /(x2 + a22/2) 



{x^ + y"^ ) d x^ ^{x^ + a? y"^) dy^ 



a;2-2Ma; + a2 + g2 + + 



Spacetime: Kerr-Newman (x, |/, ^, (KN-Euc2) 



Input file 



Contravariant tetrad 



npupkn2 . mpl 



;a;2+a2)V2 



\/2a 



2 Va;2-2Ma; + a2 + g2 ^2.2 + ^2^/2 ' 2 ^a;2 + a2 2/2 Va;2 - 2 M a; + a2 + g2 ' 
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( + a2 ) V2 



,0 



%/2c 



.2 Va:2 -2Mx + a2 + Q2 ^^.2 4, ^2 ^2 '2 + y2 ^a;2 - 2 M a; + a2 + Q2 ' 
1 ^/2 v/a;2 - 2 A'f a; + a2 + g2 



TO 



•\/a;2 + a2 j/2 



V2 



1 V2av/r^ 1 

2 ■\/a;2 + a2 y2 2 ■y/a:2 + a2 y2 -^1 — ^2 

1 \/2av/l-2;^ 1 \/2 



,0, 



,0, 



1 N/2v/-l + y2 

2 ■\/a:2 + a2 y2 
1 \/2v/-l + 2/'" 



2 ■\/a;2 + a2 ?/2 2 ^ x'^ + y"^ 1 — y"^ 2 ^ + j/2 
Corresponding line element: 
,2 (a2y2+a;2_2Mx + Q2) rft2 a ( -2 Mx + Q2 ) ( _i + ) ( ^ + 1 ) ^ ^ 

us = — h 2 



a;^ + a;-^ + a-^ t/^ 

+ ( x"* + a:2 + 2 a2 M a; - a2 Q2 _^ ^2 ^2 ^2 _ 2 y2 ^ ^ 4_ ^4 ^2 j^2 g2 ^ 



(-l + y)(y + l) rfz^^ /(x-2+a2j^2) 



( X + o y ) dx 



{x^ + a y ) dy 



a;2-2Mx + a2 + Q2 + 



Input file 



Covariant tetrad 



npdnkn2 . mpl la — 



nia 



1 \/2v/a;2-2Afx + a2 + g2 

2 ■y/a;2 + a2 y2 

1 \/2 v^a;^ + a2 y2 ^" 

2 ^a;2 -2Afx + a2 + g2 ' 

1 \/a;2 - 2 Af X + a2 + g2 

2 \J x^ ^ y2 

1 \/2 v/a:^^ + a2 y2 

2 ^x2 - 2 Af X + a2 + . 



1 V2 a v/a;2 - 2 Af x + a2 + g2 ( 1 - y2 ) 

2 ■\/x2 + a2 2/2 



1 %/2 a ^x2 - 2 Af X + a2 + g2 ( 1 - y2 ) 

2 •\/x2 + a2 2/2 



,0 



1 V2a 1 V2v/r"^(x^ +a^) ^ 1 \/2 Vx2 + y^ 

2 ■\/x2 + a2 y2 2 ■\/x2 + a2 y2 2 -^Z— 1 + ?;2 



1 V2av/r^ 1 x/2v/r^(x2+a2) ^ ^1^2^^^^ 

2 •\/x2 + a2 2/2 2 Y^x2+a22/2 2 y— 1 + j/^ 
Corresponding line element: 

^^2 _ (a^y^+a^^-2A//x + g2) dt^ ^ ^ a ( -2 A^ x + g2 ) ( -1 + ?/ ) ( y + 1 ) i z 



x2 + a2 J/2 x-^ + a'^ y^ 

+ ( x'' + x2 + 2 a2 A/ X - a2 g2 + y^ x^ - 2 Af x + a" y2 + y2 g2 ) 

[ i + y)[y + L)az /[x+ay) ^2 ^ 2 M x + + ^ (-1 + y) {y + 1 



Spacetime: Kerr- Newman (Boyer-Lindquist coordinates) (KN-BLl) jl8 



Input file 



Contravariant tetrad 



npupknS.mpl I"' — 



,1, ,0- 



r2 - 2A^r + a2 + g2 ' '' r2 - 2 Af r + a2 + g2 

1 1 r2 - 2 Af r + + g2 ^ 1 

2 r2 + a2 cos(6')2 ~ 2 



0- 



r2 + a2 cos(6')2 2 r2 + a2 cos( 6* ) 
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1 {Ism{9)r + sm{e)acos{e))aV2 1 {r ~ I acos{e)) ^2 

2 r2 + a2 cos(6l)2 ' '2 r^+a^cos{0)^ ' 
1 /(r-/acos(6'))\/2 



2 sin(6i) (r2 + a2 003(6* )2 

1 (-/sin(6')r + sin(6')acos(6'))aV2 1 (r + /acos(6') ) V2 

2 r2 + a2 cos(6')2 ' '2 r2 + a2 cos( 6* )2 ' 
1 /(r + /acos(6l)) V2 



2 sin(6') (r2 +a2cos(6')2) 
Corresponding line element: 



_ (a2 cos( 61 )2+r2 + Q2_ 2 Mr) ( -1 + cos( ))( cos( 6i ) + 1 )( Q2 _ 2 Mr ) a dt d 



r2 + a2 cos( 6 )2 
(r2 + a2cos(6l)2) d r^ 



r2 + a2 cos( y 



'cos{ef) de^ +(-l + cos(6l)) 



r2-2Mr + a2 + Q2 V 
{ cos( 61 ) + 1 ) 2 a2 cos( 9 )2 M r + a* cos( 61 )2 + r2 cos( 6* )2 
+ a2 g2 cos( 61 )2 - a2 q2 ^ 2 a2 M r + + a^^ d<^^ j {f + cos( 61 )2) 



Input file 



Covariant tetrad 



npdnknS . mp 1 I a. 



r2 + a2 003(6* )2 
'~r2-2A/r + a2 + Q2 

1 r2 - 2 M r + a2 + Q2 y 

2 r2 + a2 cos(6l)2 



,0, -asin(6')2 
^ 1 a(r2 -2Mr + a2 + Q2)sin(gl)2 

'2' '^2 r2 + a2 cos(6')2 



1 (/sin(g)r + sin(g)acos(g))aV2 ^ 1 / . u /9 

2 r2 + a2cos(^)2 , 0, - - ( r - / a cos( ) ) ^2, 

1 ( / sin( 61 ) r + sin( 9 ) a cos( 6* ) ) ( r2 + a2 ) V2' 

2 r2 + 02 003(6* )2 

1 (-/sin(g)r + sin(g)aoos(g))aV2 ^ 1 ..^79 

2 r2 + a2oos(6')2 ' ' 2 ^ v ; ; v 
1 (-/sin(6')r + sin(6i)aoos(6i)) (r2 +a2)\/2' 

' 2 



r2 + a2 oos( 6* )2 
Corresponding line element: 



_ (a2 cos( 61 )2 +r2 +g2 - 2 Mr) di^ ( -1 + cos( ))( cos( 6* ) + 1 )( Q2 _ 2 Mr ) a dt d 



(r2 + a2cos(6l)2) dr^ 
- 2 M r + a2 + g2 



r2 + a2 cos( 9 )2 



r2 + a2 oos( 9 )2 



+ (- 



-r — a cos 



(6')2) d6l2 +(_i + cos(6l)) 



( oos( 61 ) + 1 ) 2 a2 oos( 9 f M r + a* cos( )2 + r2 cos( 9 )2 

+ a2 g2 cos( 61 )2 - a2 g2 + 2 a2 M r + + a^^ d <j)^ / {r^ + cos( 9 )2) 



Spacetime: Kerr- Newman (Eddington-Finklestein coordinates) (KN-EFl) IjlQ 



Input file 



Contravariant tetrad 



npupkn4.mpl P = [0, 1, 0, 



,0, 



r2 + a2 cos( 9 )2 
IV2 



+ Ir^ — 2mr + a+C 

r2+a2cos(0)2 2 r2 + a2 cos( )2 

1 /x/2asin(6') 1 x/2 1 

2 r + /acos(6') ' ' 2 r + /aoos(6') '2 ( r + / acos( 6I ) ) sin( 6* ) 
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1 I^/2aslT^{6) 1 



1 



Iy/2 



2r-Iacos{0)' ' 2 r - I acos{0) ' 2 {r - I acos{e))sm{9) 
Corresponding line element: 



ds 



[a?coa{6Y ^-r"^ -2mr + Q'^) d-i 



2 du d r 



+ a? cos( 9 Y 

a ( -1 + cos( e)){ cos( 6) + l){~2mr + Q'^) du d (j) 



+ cos( 9 Y 

+ 2 (-1 + cos( 61 ))(cos( 61 ) + l) a dr d cf) + {-r"^ - cos{ef) dO^ + 
(-1 + 008(6*)) (cos(6l) + 1) (r^ a2cos(6')2 - 2 m r cos( 6* )2 + cos( 6* )2 

+ a^Q'^cos{9f + 2a^mr + r^ -a^Q"^ +r^a^^ d <j)^ j {r^ + cos{9f) 



Input file 



Covariant tetrad 



npdnkn4.mpl la = 



1, 0, 0,(-l + cos(6i))(cos(6') + l)a] 

1 - 2TOr + + 

2 r2 + a2 cos(6')2 ' ' ^' 

1 a ( -1 + cos( 9)){ cos( 6* ) + 1 ) ( r2 - 2 m r + a2 + Q2 y 

2 r2 +a2 cos(6l)2 
I lV2asm{9) ^1, ^ ^ 1 /V2sin(6')(r2 + a2) 

2 . + /acos(^) 'Q'2(-- + ^"-"-(^))^''2 . + 

1 ^^-^"(f\,0,-l(. + /acos(.))72,-l^^-"^(^n^-+^- 

2 -r + /acos(0) ' '2^ ^ ^ ^ ' 2 -r + /acos(0) 
Corresponding line element: 

2 (a2cos(6')2 + r2 - 2TOr + Q2) du^ 

ds = 5 ^ ——5 V2 du d r 

+ cos( 9 y 

2 a ( -1 + cos( 6* ) ) ( cos( 9) + l){-2mr + Q^) du d (j) 
r2 + a2 cos( 9 )2 

+ 2(-l + cos(6l))(cos(6i) + l)a dr rf </) + (-r2 - a2 cos( 61 )2) dd^ + 
( -1 + cos( 6* ) ) ( cos( 6* ) + 1 ) (^r2 a2 cos( 6* )2 - 2 a2 m r cos( 6* )2 + cos( 9 f 
+ a2g2cos(6')2 + 2a2TOr + r4-a2g2^^2^2^ ^^2 j [r"^ cos{9f) 

Spacetime: Kerr- Newman (Boyer-Lindquist coordinates, u = a cos 9) (KN-BL2) [|l8) 



Input file 



Contravariant tetrad 



npupknS.mpl 1°' = 



r"^ -2Mr + a? + Q'' 



,1,0, 



-2Mr + a2+Q2 



1 r2-2Afr + a2 + g2 1 
9 , 



2 r2 + ^2 



1 r2 + 

2 r2 + ^2 ' 2 

1 (/ x/a2 - u2 r + Va^ - u2 u) V2 1 Va^ - 1*2 ( r - /u ) %/2 

1 /(r-/u)a^/2 



2 \J d? — u^ ( r2 + u2 ) 
1 
2 



1 (-/Va2 -u2 



r + \/a2 — u2 u] 



r2 + u2 
1 /a(r + /u)\/2 



1 Va2-w2(r + /M)%/2 
' "2 r2 +u2 ' 



2 \/a? — ( r2 + u2 ) 
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Corresponding line element: 



,2 (u^ -2Mr + r'^ +Q^) dt^ (Q^~2Mr){a~u){a + u)dtd 
as = 7\ 2 



( + ) dr^ 



{I r + u) {I r ~ u) du^ 



r2 -2Mr + a2 + Q2 {a-u){a + u) 

{2%^ M r - - r"^ - Q"^ - r'^ - r"^ a? + Q"^ - 20^ M r) {a - u) 
{a + u) d(j)^ /{{r'^ +u^)a^) 



Input file 



Covariant tetrad 



npdnknS . mpl la ~ 



ma 



J,2 _|_ y2 



-2Afr + a2 + Q2 



,0, - 



1 r2-2Mr + a2 + Q2 1 1 (r^ - 2Mr + + Q2 ) (^2 _ ^2 ) 

2 r2 +"^ '2' '~2 a ( r2 + u2 ) 
1 Va2 - u2 ( / r + M ) ^/2 1 ( r2 + u2 ) ^ 

7 ^ ' 77 



2 r2 + u2 2 V q2--»2 (^_^ J^) 

1 {I r'^ + ur"^ + I r + u a? ) \J - v? \/2 
' 2 ^ a ( r2 + u2 ) 

1 \^2~^(/r-M)\/2 ^ 1 (r2+u^)\/2 
"2 r2 +u2 ' 2 V^ 2T:T?(-r + /M)' 

1 {I — ur^ + I r " ua"^) \J a?' — \/2 



2 a{r^+u'^) 
Corresponding line element: 
,2 - 2Mr + r2 +Q2) ^ ( Q2 _ 2 Mr ) ( a - u ) ( a + u ) rf f rf 

as = 2 



r2 + u2 a ( r2 + u2 ) 

{r^ + u'^) dr^ [I r + u) [I r — u) du^ 



r2-2Mr + a2 + (52 [a - u){a + u) 

{2u'^ Mr - a^u^ - r"^ ~ Q'^u^ - r'^ - r"^ + Q'^ -20^ Mr) {a - u) 
{a + u) d(j)^ /{{r"^ +u^)a^) 



Spacetime: Kerr- Newman (Eddington-Finklestein coordinates) (KN-EF2) IjlQ 



Input file 



Contravariant tetrad 



npupkne.mpl r= [0,1,0,0] 



1 r2 



2 mr + a2 + (52 



r2 + a2 cos( 6* )2 
/\/2 1 



,0, 



i2 cos(6l)2 
V2 



r2 +a2cos(6')2 ' 2 

1 /V2asin(6') 1 

2 r + /acos(6l) ' '2 /r-acos(6i) '~ 2 ( / r - acos( 61 ) ) sin( 6* ) 

1 /V2asin(6l) 1/^2 1 V2 



2 r-/acos(6i) ' ' 2 /r + acos(6') ' 2 ( / r + acos( ) ) sin( 6* ) 
Corresponding line element: 



2 (a2cos(6')^ +r2 - 2TOr + (52) 
rfs = 



2 du d 



r2 + a? cos( 6* )2 
2 a ( ~1 + cos( e)){ cos( ) + 1 ) ( -2 m r + g2 ) d u d 

r2 + a2 cos( 6* )2 

+ 2(-l + cos(6'))(cos(6i) + l)a dr d (j) ^ {-r'^ - a"^ cos{6f) 
(-l + cos(e')) (cos(6l) + 1) (^r2 a2cos(6')2 -2a2mrcos(e')2 + a4cos(6')2 
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Input file 



Covariant tetrad 



npdnkn6 . mpl la 



1, 0, 0, (-l + cos(0)) (008(6*) + l)a] 



1 r2 
2 



2 m r + a + I 



,1,0, 



r2 + a2 cos(6')2 

1 o ( - 1 + cos( 6* ) ) ( cos( 6* ) + 1 ) ( - 2 TO r + a2 

2 r2 + a2cos(6')2 



0') 



1 /V2asin(e') ^ 1 

2 r + /acos(6') ' ' 2 



/acos(0))\/2, 



1 V2sin(6l)(r^+a^) 

2 /r - acos(6') 



1 /V2 a sin 61 ,1, , , rr^ 1 V2 sm 6* + 
—r-T , , ( r + / a cos( 6* ) ) V2 , - -^^ , , ' 

2 -r + Iacos{e)' '2^ ^ " '2 -/r-acos(6l) 

Corresponding line element: 

{a} cos( )^ + — 2 TO r 



ds 



+2 



'2) du^ 

- 2 < 2 I a\2 V2 du d r 

+ cos( y 

a ( -1 + cos( 6* ) ) ( cos( 6') + l)(-2TOr + Q2)c;ud(^ 



+ cos( 9 Y 

+ 2(-l + cos(6l))(cos(6') + l)a dr + (-r^ - cos( 6* )2) dO^ + 
( -1 + cos( 6* ) ) ( cos( e) + l){r'^a^ cos( 6* )2 - 2 to r cos( Of cos( 6* )2 

+ a2g2cos(6')2 + 2a2TOr + r4-a2g2+r2a2) d 4>^ / {r^ + cos{ 9 )^) 



C.3 Mixmaster spacetime (Table 3) 

For this set of tests, not only are the basis vectors varied, but also their inner product. The input used for 
these tests is listed below. 



Spacetime: Mixmaster (Mix) | 
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Input file Frame 



mix. mpl Inner product of basis vectors: 

^(a)(6) ^ jjg^g^g(2a(t))^e(2b(t))^e(2c(t))^_e(2a(t)+2b(t)+2c(t))) 

Basis vectors: 

LUia = [ cos( -0 ) , sin( 2p ) sin( 9) ,0,0] 
<jj2a = [ sin( ip ) , -cos( ip ) sin( 9) ,0,0] 

W3a = [0, COS(0),1,O] 
W4a = [0, 0, 0, 1] 

Corresponding line element: 

ds2= (cos(^)2(e^(*))2 + sin(V')'(e^(*))2) rf^^ 

-2sin('(/')cos(V')sin(6l)(e'^(*)-e^(*))(e^(*)+e^(*)) d9 d (j) + 

(sin( V f sin( 6» )2 ( e^( * ) )2 + cos( V f sin( 6» )2 ( e^^ * ) )2 + cos( 6» )2 ( e'^( * ) )2) ^ 

^2 +2(e'=(*))2cos(6i) V +(e'^(*))2 di/-^ - (e^(*) )2 (eb(t) )2 (ec(t) )2 ^t^ 



Spacetime: Mixmaster (Mixl) 
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Input file Frame 



mixl .mpl Inner product of basis vectors: 

^(a)(6) ^ jig^g^g(2a(T))^g(2b(T))^g(2c(T))^_g(2a(T)+2b(T)+2c(T))^ 

Basis vectors: 



W2a 
W3q 





^2 




02 







02,0,0 



0,9. 



1 



,-Vi-e2Vi-«'2,o,o 



,0 



W4a = [0, 0, 0, 1] 

Corresponding line element: 

2_ ((e-(I^))2_(e^(^))2$2^*2(gb(T))2) ^02 



(is 



(e-i)(e + i) 

+ 2*^-(1'-l)(* + l)(eb('^) - e^(^))(e'^(^) + e^('^)) d $ + ( 

(ea(^) )2^2 _ (ea(T) )2^2 02 _^ ( gb( T ) )2 _ ^2 gb( T ) ^2 _ (-gbCT) -)2 02 

e(e^(^))2 d$ d * 



(eMT))202^2^02(gC(T))2 ^3 



(e^(^))2 rf^g 

(vE'-i)(* + i: 



(ga(T))2(gb(T))2(gc(T))2 ^ 



Spacetime: Mixmaster (Mix2) [pO|] 
Input file Frame 



mix2.mpl Inner product of basis vectors: 

jyWW =diag(l, 1,1,-1) 
Basis vectors: 

= [e^(*)cos(V), e^(*)sin(i/;)sin(6l) ,0,0] 
cj2a = [e'=(*)sin(i/'),-e''(*)cos(V')sin(6'),0,0] 

= [0, e'=(*)cos(6') ,e'=(*) ,0] 
c^4a= [0, 0, 0, e('^(*)+^(*)+^(*))] 
Corresponding line element: 

(cos(V)'(e'^(*))2+sin(^)2(e'^(*))2) 

- 2sin(V')cos(V')sin(6') (ei^^*' _ e'^C*) ) (e'^f*) + e'^^*) ) dO d 4> + 
(sin( V )2 sin( 6» )2 ( e'^( * ) )2 + cos( V )^ sin( 6* )2 ( eb( * ) )2 cos( 6* )2 ( e'=( * ) )2) d 
+2(e=(*^)2cos(e') d(l) d 4, +(e'=(*))2 _ ( g^l O )2 ( t ) )2 ( e=( t ) )2 



Spacetime: Mixmaster (Mix3) 
Input file Frame 



mix2.mpl Inner product of basis vectors: 

jyWW =diag(l, 1,1,-1) 
Basis vectors: 



VT~e2 



,e'^(^)^'Vi-e2,o,o 
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W2a 



,_eb(T)^l_02^1_^2 0,0 



,0 



a;4a = [0,0,0, e(-(^)+''(^)+^(^))] 
Corresponding line element: 

(_(ea(r) )2 + (e^(^) )2^2 _ ^2 (gb(T) )2) ^02 



ds 



(e-i)(e + i) 

2 4'^-(1'-l)(* + l)(e^ 



a(T) 



(ea(T) )2^2 _ (e^(^) )2*2 02 + (eb(r) )2 _ ^2 ^gb(T) ^2 _ (-gMT) ^2 02 
+ (eb(T))2 02^2^02(gc(T))2A ^^2 Qj£^^}^d ^ d * 



(ga(T))2(gb(T))2(gc(T))2 ^ ^.2 
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D Sample input file npupknS.mpl and sample of output 
knel .ms 



The following is the input file: 



Ndim_ 

xl_ 

x2_ 

x3_ 

x4_ 

etal2. 

eta34. 

bll_ 

bl2_ 

bl4_ 

b21_ 

b22_ 

b24_ 

b31_ 

b33_ 

b34_ 

b41_ 

b43_ 

b44_ 

Inf o_; 



4 

t 
r 
u 

phi 



(r-2+a-2)/(r-2-2*M*r+a-2+q-2) : 
1 : 

a/(r'2-2*M*r+a'2+Q~2) : 
l/2*(r-2+a-2)/(r-2+u-2) : 
-l/2*(r'2-2*M*r+a~2+q~2)/(r~2+u~2) : 
l/2*a/(r-2+u-2) : 

l/2*(I*(a-2-u-2)-(l/2)*r+(a-2-u-2)-(l/2)*u)*2-(l/2)/(r-2+u-2) 
l/2*(-(a-2-u-2)-(l/2))*(r-I*u)*2-(l/2)/(r-2+u-2) : 
l/2*I*(r-I*u)*a/(a-2-u-2)-(l/2)*2-(l/2)/(r-2+u-2) : 
l/2*(-I*(a-2-u-2)-(l/2)*r+(a-2-u-2)-(l/2)*u)*2-(l/2)/(r-2+u-2) 
l/2*(-(a-2-u-2)-(l/2))*(r+I*u)*2-(l/2)/(r-2+u-2) : 
-l/2*I*a*(r+I*u)/(a-2-u-2)-(l/2)*2-(l/2)/(r-2+u-2) : 
= ' Contravariant NPtetrad for Kerr-Newman metric 
(u=a*cos (theta) to Boyer-Lindquist coordinates)': 



The following is the annotated input/output of knel .ms. 



> restart: 

> readlib(grii) : 

> grtensorO; 



GRTensorll Version 1.26a 
December 1 , 1995 
Developed by Peter Musgrave, Denis Pollney and Kayll Lake 
Copyright 1994 ^ 1995 by the authors. 



Latest version available from : http : // astro. queensu.ca/ ~ GRHome.html 
To initiate help type Igrtensor 

Defaults read from C : \MAPLEV3\lib/ grtensor.ini 
We now load the contravariant tetrad, 
calculate the scalars, and factor them. 

> qload(npupkn5) ; 

Default spacetime = npupknS 
For the npupkn5 spacetime : 
Coordinates 
x ° = [tru(j)] 
Basis inner product 
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V 



(a) (&) - 



n = 



10 
10 
-1 
0-1 

Null tetrad {contravariant components) 

10- 



r2 - 2 M r + a2 + Q2 j.2 _ 2 M r + + Q"^ 
Ir^+a^ 1 r2-2Mr + a2 + 02 1 ^ 



2 r2 + w2 



r2 + u2 



0- 



2 r2 + u2 



1 {iVd 



2 — 1(2 r + va 



[2 r2 + u2 

1 I{r-Iu)aV2 



— {r — lu) \/2 
~ 2 ^^2— ^l2 



2 \/a^ — V? ( r2 + ^2 ) 



mbar"' = 



1 (— / Va^ — u2 7- -|- ^/a^ 



2 r2 + u2 

1 /a(r + 7^i)^/2 



_V2 1 y/a^ -u^{r + Iu)V2 
2 r2+u2 



2 \/(l2 — ^2 ( 7^2 ^ y2 y 

Contravariant NPtetrad for Kerr — Newman metric {u = a* cos{theta) to Boye\ 
r — Lindquist coordinates (J. Math. Phys. 8 265)) 



grcalc(RicciSc,WeylSc) ; 
gralter(_, radical, factor) ; 
grdisplay(_) ; 



For the npupknS spacetime : 
$00 = 
$01 = 
$02 = 



$11 = - 



*2 



2 (r + /u)2(r-/M)2 
$12 = 

$22 = 

A = 

*0 = 

*1 = 

Q^- Mr- I u M 
(r + I u) {r ~ I u)^ 
*3 = 
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*4 = 



We now generate the metric and simplify it. 

> grcalc (g(dn,dn) ) ; 

> gralter(_, simplify, factor) ; 

> grcalc (ds) ; 

> grdisplay(_) ; 



The coordinate components of Ricci and Weyl are calculated and simplified. 
(Display deleted for this appendix.) 

> grcalc(g(up,up)) ; 

> gralter(_, radical, expand, factor) ; 

> grcalc (R(dii,dn) ,Ricciscalar) ; 

> gralter(_, expand, f actor) ; 

> grcalc (C(dn,dn,dn,dn) ) ; 

> gralter(_, expand, factor) ; 

The following lines define the Weyl Scalars. 

> grdefine('P0',O,'-C{a b c d}*NPH-a}*NPm{-b}*NPl{-c}*NPm{-d}' ) ; 
Created definition for PO 

> grdefineCPl' ,{}, '-C{a b c d}*NPl{-a}*NPn{-b}*NPl{-c}*NPm{-d}' ) ; 
Created definition for PI 



> grdefine('P2' ,{}, '-C{a b c d}* (NPl{-a}*KPn{-b}*NPl{-c}*NPn{-d}-NPl{-a}*NPn{-b}*NPm{-c}*NPmbar{-d})/2' ) ; 
Created definition for P2 



> grdefine('P3',0,'-C{a b c d}*NPn{-a}*NPl{-b}*NPn{-c}*NPmbar{-d}') ; 
Created definition for P3 



> grdefine('P4' ,0, '-C{a b c d}*NPn{-a}*NPmbar{-b}*NPn{-c}*NPmbar{-d}') ; 
Created definition for P4 

The following lines define the Ricci Scalars. 

> grdefineCPOO' ,{}, 'R{a b}*NPl{-a}*NPl{-b}/2') ; 
Created definition for POO 

> grdefine('P01',O,'R{a b}*NPH-a}*NPm{-b}/2') ; 
Created definition for POl 



> grdefine('P02' ,{}, 'R{a b}*NPm{-a}*NPm{-b}/2') ; 
Created definition for P02 



> grdef ineCPll' ,{}, 'R{a b}*(NPl{-a}*NPn{-b}+NPm{-a}*NPmbar{-b})/4') ; 
Created definition for Pll 



For the npupknd spacetime : 



Line element 




r2-2Mr + a2 + Q2 {a-u){a + u) 



26 



> grdef ine('P12' ,{}, 'R{a b}*NPn{-a}*NPm{-b}/2 ' ) ; 
Created definition for P12 



> grdef ine('P22' ,{}, 'R{a b}*NPn{-a}*NPn{-b}/2 ' ) ; 
Created definition for P22 



We calculate and simplify the Weyl scalars. 

> grcalc(P0,Pl,P2,P3,P4) ; 

> gralter (_, factor) ; 

> grdisplay(_) ; 

For the npupknS spacetime 
PO = 
PI =0 
Q^-Mr-IuM 



P2 = 



(r + lu) {r — lu)^ 
PS = 



P4 ^0 

We calculate and simplify the Ricci scalars. 

> grcalc(P00,P01,P02,Pll,P12,P22) ; 

> gralter (_, factor) ; 

> grdisplay(_) ; 

For the npupknS spacetime 
POO = 
POl = 
P02 = 
1 

Pll 



2 (r2 +-u2)2 
P12 = 



P22 = 

Unfortunately, the output obtained from each of these methods are not in exactly the same form. Since 
P2 is smaller as measured in Maple words, this must be considered the fully simplified form. In fact, by 
performing the MapleV operation factor on the denominators of \E'2 (an operation requiring less than .1 
CPU seconds), it can be reduced to the required fully simplified form. In all tests listed in Appendix ^ 
the final results of calculation are either presented in exactly equivalent forms, or equivalent within some 
simplification operation requiring a negligible amount of time. 
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